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Abstract 

Working with nuclear magnetic resonance (NMR) in quadrupolar spin systems, in this paper we 
transfer the concept of atomic coherent state to the nuclear spin context, where it is referred to as 
pseudo- nuclear spin coherent state (pseudo-NSCS). Experimentally, we discuss the initialization of 
the pseudo-NSCSs and also their quantum control, implemented by polar and azimuthal rotations. 
Theoretically, we compute the geometric phases acquired by an initial pseudo-NSCS on undergoing 
three distinct cyclic evolutions: i) the free evolution of the NMR quadrupolar system and, by 
analogy with the evolution of the NMR quadrupolar system, that of ii) single-mode and Hi) two- 
mode Bose-Einstein Condensate like system. By means of these analogies, we derive, through spin 
angular momentum operators, results equivalent to those presented in the literature for orbital 
angular momentum operators. The pseudo-NSCS description is a starting point to introduce the 
spin squeezed state and quantum metrology into nuclear spin systems of liquid crystal or solid 
matter. 



* Electronic address: rauccais@cbpf.br 



1 



I. INTRODUCTION 



In a seminal work, R. J. Glauber 111] proposed the concept of a coherent state as an 
eiffenstate of the annihilation operator of the harmonic oscillator. This concept was extended 

nn 

to atoms as the well-known atomic coherent states (ACS) [2|, |3|]. These coherent states, used 
in many physical processes in the fields of quantum electrodynamics [4], trapped ions ^ 
and Bose-Einstein condensates (BECs) [Q], are of fundamental importance because they 
minimize the Heisenberg uncertainty [2'|. The minimization of uncertainty tends to amplify 
the quantum effects on the macroscopic scale. One physical system that manifests quantum 
eflFects at a macroscopic level is the BEC. In this connection, an interesting application of 
the ACS was developed by Chen et al. [7], who computed geometric phases of a coupled 
two-mode BEC system for an adiabatic and cyclic time evolution. Other applications are 
the generation of macroscopic quantum superpositions jsl] and dynamic control jol, [lo| in 
BECs. 

The topic of geometric phase has recently received considerable attention and it has been 



employed increasingly in quantum information processing 
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iJ]. The quantum version 



of geometric phases was first reported and studied 



3yM.V. Berry in the 1980's [15j and 



implemented, thereafter, in many physical systems jl6l-ll8|. A general discussion of geometric 
phases acquired by pure states was published by N. Mukunda and R. Simon [19]. In this new 
approach, the authors took advantage of quantum kinematic concepts, so that this formalism 
could be applied in a general context, such as in any smooth open path defined by unit vectors 
in a Hilbert space. Many applications using this formalism have appeared, for instance in 



noncyclic geometric quantum computation [20|] and, more recently, nonadiabatic geometric 
quantum gates using composite pulses [isl , both of these applications being discussed in the 
NMR context. Other applications of the Mukunda-Simon approach have been the study of 
geometric phases in a nonstationary superposition of atomic states induced by an engineered 



reservoir 



2l[ | , and also the study of the eflFect of geometric phases on the pseudospin dynamics 



of two coupled BECs [22] and of a spin-orbit-coupled BEC [23|]. Likewise, NMR experiments 
in systems with two coupled spins and quadrupolar spin systems have been published, in 
which the geometric phase of an adiabatic cyclic evolution of an initial quantum state is 



ana. 



yzed 



24l |. Other potential NMR applications of adiabatic jllL |25| and nonadiabatic 



26|-l28| geometric quantum computation have also been reported. 
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These advances inspired us to study two subjects related to NMR quadrupolar systems. 
The first is the possibility of transferring the definition of ACS to a nuclear spin state, which 
will be referred to as a pseudo- nuclear spin coherent state (pseudo-NSCS). The second 
concerns an application of pseudo-NSCSs to compute geometric phases in three diflFerent 
configurations: i) the free evolution of the NMR quadrupolar system [29! and, by analogy 
of the evolutions of the NMR quadrupolar system, that of ii) single- and in) two-mode 
BEC-like system [30, IMI]. Under this cyclic evolution, we apply the Mukunda and Simon 
formalism and obtained theoretical results analogous to those of Zhang et al. in Ref. |32| . 
Special attention was given to configurations ii) and in), in which we took the advantage 
of the mapping between the quasiparticle description and pseudo-spin momentum operators 
for any spin value /. 

In general, the present paper presents experimental results and theoretical applications 
of the mapping between the atomic and nuclear spin scenarios. This will be a basis for 
future experimental implementations of quantum control, and also to establish an NMR 
quadrupolar system as a workbench for a BEC system. 

This paper is organized as follows: In section [III the main concepts and a review of ACS 
are presented, while in section [iTT] the NMR quadrupolar system is briefiy described. Section 
[TVl is devoted to the experimental procedure to initialize the pseudo-NSCSs and describes 
the quantum control of polar and azimuthal rotations. In section [Vl theoretical applications 
of the pseudo-NSCS based on the geometric phase concept are elaborated for configurations 
i), ii) and Hi) and in section IVTl we report our conclusions. 



II. ATOMIC COHERENT STATES 



In the theory of ACS, a Bloch vector n = (sin 9 cos 99, sin 9 sin (f^ cos 9) is transformed by 
the action of a rotation operator R^i^^ = e~^^^'^'^\ The operator R^i^^ represents a rotation of 
angle 9 around an axis defined by m = (sin 99, — cos Lp^ 0), where J • m = J^^ sin (/^ — cos (f. 
The operator J signifies the total angular momentum, with components J^;, J^, and J^, such 
that j is the total angular momentum quantum number and its projection corresponds to 
m = — j, — j + 1, . . . , j — 1, j. After a simple algebraic procedure, the application of He^cp to 
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FIG. 1: (Color online) For a nuclear spin system we define rotation angles 9 and Lp on the Bloch 
sphere relative to the xyz frame and the corresponding value of Q in the complex plane denoted 
by the Greek capital form of The north pole corresponds to ^ = and the south pole to 
C = oc. 



the fundamental state denoted by Ij, —j) produces an excited state represented by 

1 



-e 



(1 + cxy 



li, -j) , 



j-\-m 



(i + C*C) V (i + "^)! (i - "^)! 



(1) 



where C = — e^^'^tanl, with the angles if and 9 as shown in Fig. 1, while J± = J^^ ± zJ^, 
and Ij, m) are eigenstates of the operator J^, with eigenvalues mh. The complex phase that 
characterizes the ACS may be defined in the complex plane Z (Greek capital Q or over the 
Bloch sphere, where the ground state is represented by the south pole and the most excited 
state is identified with the north pole. 

Let us consider the density operator of the ACS denoted by p = \C{9^ip)) {C{9^ip)\ = 
Y.in,m'=-j Pm,m' |i, ^) 0', rn'\. From Eq. ([I]), each element of p is given by 



^^2j+mWgiK-m)^ ^^^2j-m-m' (^/2) sin2i+^+^' [9/2) 



[j + m') \ [j — m')\ y [j + m) \ [j — m)\ 

Hence, if the angular parameters 9 and Lp are known, the amplitude of each element 
Pm,m' can be obtained. The inverse procedure, that is, the computation of 9 and Lp from the 
elements Pm,m'-> is also possible. In the case of 0, we use 
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= sm 



cos 



2' 



2' 



(3) 



for any value of j. In the case of when j is an integer 



(Po.i+Pi,o) + (])<■ ^ I,.. 

2,(-lf'«(cosf)='-'(smf)''*'(2,)! 

Similar expressions can be obtained in terms of the density matrix elements po,-i and p-i,o 
(not shown). For j a half integer, these equations become 



2(-l)2^(cosf)^^(smf)^^(2j)! 



cos(/p, (6) 



= sin(/;. (7) 



2^(-lf (cosf)^^ (sin|)'^'(2j)! 
Therefore a general ACS given by ([I]) is completely characterized by 6^, and Lp. 



III. NMR QUADRUPOLAR SYSTEMS 

NMR quadrupolar systems are composed of nuclei with spin / > 1/2 subjected to both a 
magnetic field and an electric field gradient, where the nuclear magnetic moment is quantized 
as M = I^I — 1,...,—/. Considering also the interaction with a radio frequency {RF) field 
used for excitation, in a reference frame rotating around z-axis with angular velocity ujrf^ 
the NMR Hamiltonian can be written (for more details see section VII of Ref. {29!): 

Hnmr = -h {ul - ujrf) + (31^ - I^) 

+huji (I^ cos Lps + ly sin LPs) + T-Lenv^ (8) 

where I^;, I^, are the x^y^z components of the total spin angular momentum operator. 
The first term of the Hamiltonian is due to the interaction of the nuclear magnetic moments 
with the strong magnetic field Bq (Zeeman term). The second is due to the interaction 
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of the quadrupole moments of the nuclei with the internal electric field gradient, {ul) 
and {(jUq) stand for the Larmor and Quadrupolar frequencies (coupling strengths), which in 
this case, satisfy the inequality |cl;q| <C The third term represents the external time- 

dependent RF field perturbation of intensity Bi =6^1/7, applied along the direction of the 
corresponding vector m, with 7 being the gyromagnetic ratio. The fourth term represents 
weak interactions between quadrupolar nuclei and other nuclei, electrons, random fields, 
referred to here as environment. This term will be neglected as its contribution is weak {34]. 

A quadrupolar system in the rotating frame will be used to describe two different situa- 
tions. In the first, the external time-dependent perturbation is present, with \(jJl\ = \ujrf\^ 
and \(jJi\ » \(jJq\ specifies a regime where the perturbation equally affects all energy levels 
of the system, non-selective pulse; the evolution operator is then given by 



Urf = exp [-iujit (I^ cos ips + T-y sin (fs)] • (9) 

The physical parameters to be identified are the polar (9) and azimuthal {(f) angles defining 
the rotation operator R6i,(^. From the RF operator Urf (t), these angles are 



9 = cjit. 



TT 



(10) 



Secondly, when the external RF perturbation is absent, the corresponding evolution 
operator is 



U 



Ev 



exp 



iK-c.HF)tI.-^ (31^-1^) 



(12) 



which refers to a free evolution regime. Note that the operator Uev cannot be written as 
a rotation operator R6i,(^, since it is not generated by the vectors n and m (see section ITTl) . 
This approach is similar to that discussed by Kitagawa and Ueda jssi in the study of spin 
squeezed states. 



IV. EXPERIMENTAL IMPLEMENTATION OF THE PSEUDO-NSCS 

The NMR experiments were carried out at room temperature in a VARIAN INOVA 400 
MHz spectrometer. The pseudo-NSCS was implemented with the ^^Na nuclei (/ = 3/2) 
present in sodium dodecil sulfate (SDS) in a lyotropic liquid crystal, prepared with 21.3 wt 
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% SDS, 3.7 wt % decanol, and 75 wt % D2O. The strength of the Larmor frequency and 
quadrupolar couphngs were 105.85 MHz and 15 kHz, respectively. Typical vr-pulse lengths 
of 8 /xs and recycle delays of 500 ms were used. The T2 and Ti relaxation times of the ^^Na 
nuclear spins are 2.6 ± 0.3 ms and 12.2 ± 0.2 ms, respectively. 

Room temperature thermal equilibrium NMR systems are almost maximum mixture 
states, represented by a density matrix which deviates only slightly from the normalized 
identity matrix 

(13) 

where /3 = l/ksT and Z = Tr ^e^-f^^NMR)^ partition function, T the room temper- 

ature and kB the Boltzmann constant. Using suitable spin rotations - RF pulses and free 
evolutions - and a temporal average technique (so-called strong modulated pulses - SMP 
js^), this thermal state can be transformed into a state of the form 



P~ \^--ejl + e\ij){ij\, (14) 

where lip) {ipl is called a pseudo-pure state, having the form of a pure state density matrix 
with unitary trace, and e = /3hujL/Z. By the technique of quantum state tomography 
(QST) I37I, we obtain experimentally the quantum state lip) that will be used to represent 
the pseudo-NSCS |/,/) = |C(0,0)). This fundamental pseudo-pure state represents, from 
the NMR point of view, the precession of the magnetic moment around an axis defined by 
the direction of the strong static magnetic field Bq. To obtain other excited pseudo-NSCSs, 
we just need to apply the operator Urf (Eq. ([9])) to |/,/). In Fig. [2] we show how to 
implement the fundamental and excited pseudo-NSCSs. 



A. Initializing pseudo-NSCSs 

To produce excited pseudo-NSCSs from |/,/), we apply non-selective pulses of strength 
uji and phase (fg during time interval t = O/uJi (see Eq. ([9j)). As a first example, the excited 
pseudo-NSCS |C (7r/2, 0)) = | — 1) can be generated by applying a non-selective 7r/2 pulse that 
produces a rotation about the negative ?/-axis. On the other hand, a rotation of 7r/2 about 
the positive ^-axis of the vector state |/, /) generates the state \( (7r/2, tt)) = |1). We observe 
that the unitary transformation (of Eq. ([9])) associated with non-selective pulses, acting on a 
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FIG. 2: (Color online) The pulse sequence has four steps. In the first step the pseudo-NSCS |/, /) 
is implemented by the SMP technique. The second step consists of the non-selective pulses with 
parameters {0^(p) that produce excited pseudo-NSCSs \({9^Lp)). The third step is the azimuthal 
rotation through time r, which can be used to produce cyclic evolutions. The last step involves 
the read-out of the pseudo-NSCS (QST procedure). 

given pseudo-NSCS, produces excited pseudo-NSCSs analogously to those discussed in Ref. 



In Fig. 3, we present experimental results for the deviation density matrix, in the form 
of bar charts of the real (left) and imaginary (middle) components, for the quantum states 



\aO,^)) = |0), |C(vr/2,7r)) = |1), |C(7r/2,0)) = |-1), |C (7r/2, 37r/2)) = |^), |C (7r/2, 7r/2)) = 



|— i), and \( (tt, (f)) = |oc) (see Tab. IIV Aj) . The intensity of the elements of each deviation 
density matrix was obtained by QST. On the right of Fig. 3, the Husimi Q-distribution 
was used to characterize these implemented quantum states, which is sketched on the Bloch 
sphere [sl. The highest positive normalized intensity corresponds to the red color and the 
most negative intensity corresponds to the orange color. The highest intensity (red color) 
of the quasi-probability distribution function indicates the orientation of the spin nuclei 
magnetization. The large dispersion of the Husimi Q-distribution function is due to the 
small spin value / = 3/2, or analogously, the small number of particles = 2/ = 3 of an 
atomic system. 

In order to compare the theoretical predictions with the experimental values of the de- 
viation density matrices of Fig. 3, we present in Tab. IIV Al the complex phase ( (^, ip) and 
the Bloch vector n. The experimental values for 9 and (p were computed with Eqs. (|3j), (|6]), 



and ([ZD. 
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FIG. 3: (Color online) Experimental results for the ^^Na nuclei in the sample of SDS. The QST 
results are represented by bar charts (left and middle column) for the implemented pseudo-NSCSs 
and on the Bloch Sphere (right column) by the Husimi Q-distribution. The top bar chart represents 
the pseudo-NSCS |0), the left (right) bar chart showing the real (imaginary) part of the deviation 
density operator, and similarly for the other bar charts. The pseudo-NSCS |0) is implemented as 
explained in section [IVl by the SMP technique. The bottom bar chart shows the pseudo-NSCS 
|oo), a state that is implemented by applying a rotation operator R^^^^ with 9 — n and Lp — 'iTi/2 
to the pseudo-NSCS |0). The pseudo-NSCS | — 1) is implemented by applying a 7r/2 non-selective 
pulse in the negative y— direction to transform |0). Analogous remarks apply to the bar charts 
labeled as |1), and In the right column, we present the calculated Husimi Q-distribution 

of each deviation density matrix and, at the bottom, the intensity of the Q-distribution function 
is encoded in the colour bar. 




B. Controlling pseudo-NSCSs 



Besides knowing how to initialize a quantum system, i.e., how to prepare the desired 
quantum state, an important task in quantum physics is the control of the quantum system 
with high accuracy. The implementation and control of pseudo-NSCSs are important for 
applications in quantum computing and geometric phases [22], squeezed states [35], and 
quantum metrology 381]. In the present case, such control is performed through polar and 
azimuthal rotations. 

Polar rotations - In the context of polar rotations, we have demonstrated the experimen- 
tal control of angle 9 by taking the initial quantum state |C(0,0)) = |0) and applying the 
rotation operator R6i,(^, with 9 = 0, 7r/18, 27r/18, . . . ,7r (total of nineteen steps) and fixed 
Lp = TT. After each step, the state \( (^, ip)) was reconstructed by QST (see pulse sequence in 
Fig. 2). First, the average values of the magnetic angular momenta {ix,y,z) were calculated 
for each deviation density matrix \( {9^ (f)) {( {9^ The experimental results (symbols) 
are shown in Fig. 4 together with the values obtained by numerical simulation (solid lines) 
and theoretical prediction (dashed lines). We use the term numerical simulation when we 
apply numerical recipes to the dynamics of evolution given by Hamiltonian ([8j), while the 
theoretical prediction is given by the evolution operator of Eq. ([9]) . The differences between 
simulated data and theoretical prediction are caused by the quadratic term of Hamiltonian 





Theoretical 


Experimental 


e 




IC(^,^)> 


n 


IC(^,^)) 


n 








|0) 


(0,0,1) 


|0.02 - O.OH) 


(-0.05,-0.02,0.99) 


TT 

2 





|-i> 


(1,0,0) 


1-0.88 -0.12i) 


(0.98,-0.14,0.12) 


TT 

2 


Stt 
2 




(0,-1,0) 


|-0.88i) 


(0, -0.99,0.13) 


TT 

2 


TT 


|i) 


(-1,0,0) 


|0.79 + 0.25i) 


(-0.94,0.30,0.18) 


TT 

2 


TT 

2 


\i) 


(0,1,0) 


1-0.09 + 0.84i) 


(0.1,0.98,0.16) 


TT 





oc) 


(0,0,-1) 


|7.75 + 1.16i) 


(-0.25,0.04, -0.97) 



TABLE I: Comparison between theoretical and experimental values obtained from the deviation 
density matrix of Fig. 3 by means of equations ([3]), ([6j) and ([7j). Experimental values are rounded 
to two decimal places. 
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FIG. 4: (Color online) Experimental results (symbols), numerical simulations (solid line), and 
theoretical predictions (dashed lines) of a polar rotation of the Bloch vector for the initial pseudo- 
NSCS |0), implemented experimentally by the SMP technique. QST is performed for nineteen 
different values of 9 from to tt. The average values of the magnetic angular momentum components 



(\x,y,z) and the total angular momentum (I) = \jYlii=xy z (-'-^)^ shown in a). The accuracy of 
the polar rotations is calculated by the fidelity F (Eq. (IT5]l ) between experimental results and 
theoretical prediction, as depicted in b). 



dHj). Such a term is the main source of imperfections in the generation and control of the 
pseudo-NSCSs. 

In order to find the fidelity of the experimental implementations, we compared the tomo- 
CTaphed deviation density matrix pExp with the theoretical prediction prhe^ by calculating 

Tr {pExp ' PThe} 



F = 



(15) 



^Tr {pExp ' PExp} Tr {prhe ' PThe} 
The results for F are presented in Fig. 4b and show that the experimental generation and 
control of the pseudo-NSCSs can be achieved with fidelity higher than 80 %. 

Azimuthal rotations. -The azimuthal rotation is performed around the z-axis so that it can 
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be achieved by the pulse sequence [7r/2]^ i^^'ly ~^ \^ I'A-x^ where the subscripts represent 
the phases of the non-selective pulses and the values in square brackets the rotation angles. 
This pulse sequence is referred to as composite pulse jsol. As a first step in the experiment, 
we arbitrarily prepare the excited pseudo-NSCS |C (7r/2, tt)). By applying the composite 
pulse sequence to the initial state |C(7r/2,7r)) we can obtain thirty-three distinct pseudo- 
NSCSs IC (7r/2, TT + ^r)) with angular values Lp^ = 0, 7r/16, 27r/16, 37r/16, . . . , 27r and 9 fixed 
at 7r/2. Finally, the transformed state is reconstructed by the QST procedure. 

The average values of magnetic angular momenta {ix,y,z) were computed for each devi- 
ation density matrix l^" (7r/2, tt + (fr)) (C (^/2, tt + (fr)\' The experimental results (symbols) 
are plotted in Fig. 5a together with the values of the numerical simulation (solid lines) and 
theoretical prediction (dashed lines). As in the case of polar rotation, using one non-selective 
pulse, the diflFerences between simulated data and theoretical prediction are caused by the 
quadratic term of Hamiltonian ([8]). However, in the study of azimuthal rotation we need 
four pulses: one to create the excited state \( (7r/2,7r)) and three making up the composite 
pulse sequence. These additional pulses are another source of errors in the generation of 
the state, due to the error accumulated during the calibration of each pulse and time delays 
between them, which are of order 2-3 //s, despite being small, is not truly negligible. In 
the azimuthal rotation, the extra pulses and delays take 10-12 /xs more time to implement 
than a polar rotation by the same angle. During this additional time, the contribution of 
the quadrupolar coupling is more evident, transforming the implemented pseudo-NSCS with 
a non-null degree of squeeze [35|]. This is the reason why the fidelity F in Fig. 5b is, on 
average, lower than the fidelity in the polar case (see Fig. 4b). In order to diminish these 
errors it is necessary to calibrate non-selective pulses of short duration. 



V. THEORETICAL APPLICATION OF THE PSEUDO-NSCS 
A. Geometric phase 

In this section, we use the Mukunda- Simon formalism fl9^ to compute geometric phases 
that arise from the three distinct cyclic evolutions. For a system evolving according to 
the Schrodinger equation, the geometric phase (j)G can be defined as the difference be- 
tween the total phase (pr {^) = arg (0) I?/; (r))] and the dynamic phase (poi^) = 
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FIG. 5: (Color online) Experimental results (symbols), numerical simulations (solid line), and 
theoretical predictions (dashed lines) of an azimuthal rotation of the Bloch vector for the initial 
pseudo-NSCS implemented experimentally by the SMP technique. QST is performed for thirty- 
three different values of from to 27r. The average values of the magnetic angular momentum 



components (\x,y,z) and the total angular momentum (I) = \l^i=xy z (-'-^)^ shown in a). The 



accuracy of the azimuthal rotations is calculated by the fidelity F (Eq. (|T5l) ) between experimental 
results and theoretical prediction as depicted in b). 

-^r(V'(t)|||t/'(t)),i-e., 

<i>G (r) = arg [(t/; (0) \^ (r))] + i f {i^ {t)\ ^ {t)) dt, (16) 



where r is the overall time of the evolution 122 



m 



Before starting the calculation of the geometric phases, it is worth mentioning that this 
phase, acquired by the quadrupolar NMR system, is related only to the pseudo-pure state. 



since its definition (fT6l) does not hold for general mixed states |42|. Assuming an initial 
pseudo-NSCS \({9^Lp))^ the geometric phase acquired by this quantum state is calculated 
while undergoing three distinct cyclic evolutions: i) the free evolution under the NMR 
quadrupolar Hamiltonian (Eq. (|8j)), and the evolution of such a system under the analogous 
ii) single- and in) two-mode EEC Hamiltonian. This analogy between the NMR quadrupolar 
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system and BECs was explored in Ref. [30|] through the Holstein-Primakov formahsm. 



B. Free evolution of the NMR system 

The first case to be analyzed is the free evolution of the quadrupolar system with the 
offset frequency set uorf = ujl—ujq. In this case the evolution operator, Eq. ( fT2l) . is rewritten 
as 



exp 



-lip 



-h + 1(311 -I') 



(17) 



where (f^ = (jOqI^ which from now on will vary from to 27r. Using Eq. ( fT6l) and computing 
the average values of operators and I for the pseudo-NSCS |C {0^ (f))^ the geometric phase 
can be expressed by: 



Vg (C) = 7r/ 



1 



cos^ 6* - 2 COS 6* + 



(18) 



2J 2 

This expression has a quadratic dependence on the harmonic function, because the Hamilto- 
nian in this analysis is quadratic in the nuclear spin angular momentum operator. Applying 
Eq. f|T^ to a spin 1 = 1/2 nucleus, we find 

^0G(C)=7r(l-cos^), (19) 



which corroborates the well known results reported for two-level particles in references 
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431]. Note that in this case the single harmonic function arises from the linear 



dependence of the nuclear spin angular momentum operator in the Hamiltonian, which 
allows the connection between the geometric phase and half of the solid angle on the Bloch 
sphere. On the other hand, for spin / = 3/2, the presence of the quadrupolar term makes 
the geometric phase proportional to the square of the solid angle. 

C. NMR system evolving as a single- mode BEC 



As a second example of the free evolution of the NMR quadrupolar system, we explore 
its analogy with the single-mode BEC, as reported in Ref. js^. In this description, the 
frequency urf of the evolution operator in Eq. (fT2l) satisfies the mathematical expression 



Uev {^') = exp 



-4(i?-(2.-:,i.-^^: 



(20) 
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where the angular parameter is denoted by ^p^ = uoqi. The geometric phase acquired by the 
pseudo-NSCS |C ^)) transformed by the operator (l20l) is 

7r/(2/- 1) 



(C) 



(3 + cos^)(l-cos^). 



(21) 



In this interpretation, the geometric phase displays a qua-dra-tic dependence on the har- 
monic function, which is linked to the quadratic dependence of the number operator in the 



Holstein-Primakov representation 



44l |. We note that for 7 = 1/2 the geometric phase is 



null. This is consistent with a discussion on geometric phase in a study of the quasiparti- 
cle dynamics in a BEC where one particle becomes a free particle and is not part of the 



condensate 



32|. 



D. NMR system evolving as a two- mode BEC 



In the third case, the two-mode BEC Hamiltonian in the Schwinger representation (see 
Eq. (14) of Ref. [7]) is simulated by Hamiltonian ([8]) of the NMR system. In this de- 
scription, the pseudo-angular momentum operators {3i=x,y,z) are mapped on to the nu- 
clear spin momentum operators {li=x,y,z)^ and the Hamiltonian parameters are mapped by 
ujq {ujrf — ^l)) Q G uji and (f ^ (fg- We will work in a particular config- 

uration where the scattering lengths satisfy the relation Aa : Aab : = 1.03 : 1 : 0.97 for 



modes a and b of the BEC 



22| . implyin 



resonances, as discussed in references |6 



Q = 0- This is usually achieved through Feshbach 



7. 
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4l|. We also impose the condition that the 



external laser field coupling the two modes is turned oS (G = 0). The net atomic evolution is 
a rotation around the z-axis of frequency cjq, which one in the NMR scenario is implemented 
by means of the composite pulse sequence explained in section HVl i.e. 



where the pulse sequence is read from left to right. After a cyclic evolution of state |C (0, 99)), 
the geometric phase is: 



>G (C) = -27r/(l-cos6^) 



(22) 



reproducing a similar mathematical expression in Eq. (14) of Ref. j4l| for an adiabatic 



evolution. In this case the geometric phase displays a linear dependence on the harmonic 



15 



function, which reflects the hnear dependence of the Hamiltonian in the nuclear spin angular 
momentum operator. 

As in the Section IVBl the particular case of the geometric phase for a nuclear spin 7 = 1/2 
gives: 

(C) = -7r(l-cose), 

a result that resembles the previous expression f fTOl) for a spin 7 = 1/2 particle. The negative 
sign means the opposite direction to the nuclear spin rotation. Indeed, this is due to the 
equivalence between superposition of spin |t) and spin ||) and the ACS deflnition for a two 
level-system. 



VI. CONCLUSION 



In this paper we have described how it is possible to transfer the ACS concepts to the 
nuclear quadrupolar system labeled as a pseudo-NSCS. We experimentally showed how to 
initialize and to tomograph the pseudo-NSCS. The experimental manipulation of this state 
was performed by polar and azimuthal rotations, in which we note that the greatest obstacle 
to perfect initialization and control is the quadrupolar term of the NMR Hamiltonian. We 
also studied theoretically the geometric phase generated by the evolution of the pseudo- 
NSCS in three cases: i) free evolution of the NMR system, ii) single and in) two mode 
BEC. In the first and second cases, the geometric phase is a quadratic function of cos 9 for 
7 > 1/2, which is a signature of the quadrupolar term. For 7 = 1/2, the geometric phase is 
7r(l — cos^) for free evolution of the NMR system and null for the single mode BEC case. 
In the latter case, for > 0, the null geometric phase indicates there is no condensation. 
In the third case, the geometric phase acquired by the pseudo-NSCS under the action of 
a composite pulse sequence reproduces the same result as a geometric phase generated by 
a cyclic azimuthal rotation. The present study would be still more complete if one could 
measure experimentally the geometric phase discussed in each case. This task could be 
undertaken with two coupled nuclear spin systems, one spin 7 > 1/2 and the other S > 1/2 
(by analogy to Ref. |ll]), which we will try to implement in the near future. Finally, the 
concept of pseudo-NSCS in the NMR quadrupolar scenario offers new insights into spin 
squeezed states and quantum metrology implemented in liquid crystals or solid matter. 
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